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METHOD A ND APPARATUS FOR OPTICAL MEASUREMENT 
OF CONCENTRATION AND TEMPERATURE OF LIQUIDS 


STATEMENT OF GOVERNMENT RIGHTS 
This invention was made with United States Government support 
5 under National Science Foundation Contract No. CTS-9702644. The United States 
Government may have certain rights in the invention as provided in the Contract. 


FIELD OF THE INVENTION 
The present invention relates to measurement of optical properties of 
materials, and, more particularly, to precision optical measurement of refractive 
1 0 index, concentration and temperature of materials. 


BACKGROUND OF THE INVENTION 
The precise measurement of concentration in liquids is important in 
fields such as chemical analysis and processing, diagnostics, semiconductor 
manufacturing, waste inspection, and measurement of liquid difiusion coefficients. 

1 5 Measurement of liquid temperature is also important, as is determination of the 
refractive index for both liquids and non-liquid materials. As such, a variety of 
techniques have been developed to measure concentration, temperature and/or 
refractive index. Bergman et al. developed a fiber-optic probe to measure salinity 
distribution in liquids, as described in T. L. Bergman, F.P. Incropera and W. H. 

20 Stevenson, "Miniature Fiber-Optic Refractometer for Measurement of Salinity in 

Double-Diffusive Thermohaline Systems," 56 Rev. ScL Instrum 291-96 (1985) and T. 
L. Bergman, D. R. Munoz, F. P. Incropera and R. Viskanta, "Measurement of Salinity 
Distributions in Salt-Stratified, Double-Diffusive Systems by Optical Deflectometry," 
57 Rev. ScL Instrum 2538-41 (1986). Other techniques include a planar laser- 

25 induced fluorescence technique, as described in A Lozano, S. H. Smith, M. G. 
Mungal, and R. K. Hanson, "Concentration Measurements in a Transverse Jet by 
Planar Laser-Induced Fluorescence of Acetone," 32 AIAA Journal 218-21 (1993), and 
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an invasive heat-marker method as described in V. A. Vink, G. A. Sokolov, and Yu. 
S. Fomochev, "Measurement of the Concentration of Flowing Liquid Solutions," 58 
Journal of Applied Chemistry of the USSR 357-59 (1 QRS) Interferometric techniques 
are described in T. Konishi, S. Naka, A. Ito and K. Saito, "Transient Two- 
5 Dimensional Fuel-Concentration Measurement Technique," 36 A pplied Optics 88 1 5- 
19 (1997), T. A. Wilson and W. F. Reed, "Low Cost, Interferometric Differential 
Refractometer," 61 Am. J Phvs. 1046-48 (1993) and R. J. Harris, G. T. Johnston, G. 
A. Kepple, P. C. Krok and H. Mukai, "Infrared Thermooptic Coefficient 
Measurement of Polycrystalline ZnSe, ZnS, CdTe, CaF 2 , and BaF 2 , Single Crystal 
10 KCI, and 11-20 Glas," 16 Applied Optics 436-38 (1977), and a phase-locked-loop 
ultrasonic method is described in K. Ikeda, "Ultrasonic Measurement of 
Concentration in Solutions by a Phase-Locked Loop Method," 36 Jpn. J. Appl Phyx 
3180-83 (1997). 

Refractometers are routinely used to evaluate the refractive index to 

15 determine the concentration of a liquid mixture, as described in J. E. Geake, "Linear 
Refractometers For Liquid Concentration Measurement " Chemical Engineer 305-08 
(1975). Still other techniques reported in the literature for accurately measuring the 
refractive index of solids and gases include the minimum deviation method as set 
forth in L H. Malitson, "Refractive Properties of Barium Fluoride " Journal of the 

20 Opticql Swfcty of America 628-32 (1964) and B. C. Piatt, H. W. Icenogle, J. E. 

Harvey, R Korniski and W. L. Wolfe, "Technique for Measuring the Refractive Index 
and Its Change with Temperature in the Infrared," 65 Journal of the Optical Society nf 
America 1264-66 (1975), the use of a Littrow prism as described in E. D. McAIister, 
J. J. Villa and C. D. Salzberg, "Rapid and Accurate Measurements of Refractive 

25 Index in the Infrared " 46 Journal of th e Optical Society of America 485-87 (1956) 
and A. R. Hilton and C E. Jones, "The Thermal Change in The Nondispersive 
Infrared Refractive Index of Optical Materials " 6 Applied Optics 1513-17 (1967), 
Brewster angle techniques as described in I. K. Smirnov, Y. G. Polyakov and G. N. 
Orlov, "Arrangement for Measurement of Index of Refraction and Thickness of 

30 Transparent Dielectric Films by an Optical Method," Journal of the Optical F acjpf y nf 
America 546-47 (1980), and others. 
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Examples of other techniques are set forth in S. M. Chernov, K. K. 
Zhilik and P. G. Rabzonov, "Determination of the Index Refraction of Liquids and 
Gases in Capillaries," 37 Journal of A pplied Spectroscopy (English Translation of 
Zhurnal Prikladnoi Spektroskopii) 1069-72 (1982), L. A. Danisch, "Removing Index 
5 of Refraction Constraints in the Optical Measurement of Liquid level" Fiber Optics 
qndLgser $emonX26S-79 (1992) and D. R. Lide, Ed., CRC Handbook of Chemistry 
and Physics (CRC Press, Boca Raton, 1998). All of these techniques, however, suffer 
from one or more of the following shortcomings: direct contact is required with the 
material being measured, poor resolution is noted, complicated and expensive 

1 0 components are required, systems are physically large and difficult to operate, or a 
visual, subjective analysis of the data is relied upon. 

In view of the foregoing deficiencies of currently known techniques for 
measurement of concentration and temperature in materials, including liquids, a need 
clearly exists for an apparatus and method for precision, non-contact measurement of 

15 concentration and temperature of a material, such as a liquid, where the apparatus and 
method does not require complicated and expensive components, is compact, easy to 
operate, and does not rely on visual, subjective readings of measurement data. 

SUMMARY OF THE INVENTION 
The present invention, which substantially overcomes the 

20 shortcomings of the currently known techniques, provides a method for determining 
the concentration and temperature of a transparent liquid. The method includes the 
steps of causing the liquid to be contained in a vessel having a transparent entrance 
side and a transparent exit side; causing a beam of light to impinge on the entrance 
side; and then calculating the refractive index of the liquid using SnelFs law. In the 

25 step of causing the liquid to be contained in the vessel, the vessel can have an entrance 
side and an exit side with a known angular relationship therebetween, and can be 
immersed in known surroundings. In the step of causing the beam of light to impinge 
on the entrance side of the vessel, the impingement can be at an angle 0 S with the 
normal to an outer surface of the entrance side. The beam can then pass through the 
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material, and then through the exit side, from which it exits at an angle 0 e with respect 
to the normal to an outer surface of the exit side impinged by the light beam. 

The calculation of the refractive index of the material using Snell's law 
can be based on the angles 6; and 6 e , and can be done by applying Snell's law at the 

5 interface between the surroundings and the entrance side, the interface between the 
entrance side and the material, the interface between the material and the exit side, 
and the interface between the exit side and the surroundings. 

The present invention also provides a method for determining the 
concentration of a given component (e.g., a solute) in a sample of a multi-component 

10 liquid mixture (e.g., a solution). The method includes the step of determining the 
refractive index of the sample of the multi-component liquid mixture as described 
above and then comparing the determined refractive index of the sample of the multi- 
component liquid mixture to predetermined data relating different concentrations of 
the given component of the multi-component liquid mixture to corresponding values 

15 of the refractive index of the multi-component liquid mixture. In this manner, the 
concentration of the given component of the sample of the multi-component liquid 
mixture can be determined from the refractive index of the sample of the multi- 
component liquid mixture determined as described above. 

The present invention yet further provides a method for determining a 

20 change in concentration of a given component of a sample of a multi-component 
liquid mixture from an initial concentration of the given component which 
corresponds to an initial refractive index of the sample of the multi-component liquid 
mixture, at an initial sample temperature. The method includes determining an initial 
and a subsequent refractive index of the sample of the multi-component liquid 

25 mixture in the manner described above, and then determining a change in the 
refractive index by subtracting the initial refractive index from the subsequent 
refractive index. The method further includes measuring the initial and subsequent 
temperatures of the sample of the multi-component liquid mixture and then 
determining the change in the concentration of the given component in the sample of 

30 the multi-component liquid mixture from the initial concentration according to the 
following approximate formula: 
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AC » (An - (dn/ffT)AT)(dn/dC)-\ (1) 

where: 

An is the change in the refractive index of the sample, 
AC is the change in the concentration of the given component in the sample of 
5 the multi-component liquid mixture from the initial concentration, 

dn/3T is partial derivative of index of refraction, n, of the multi-component 
liquid mixture, with respect to temperature, determined from known data (can be 
evaluated at sample temperature and is approximately constant for a AT up to 10- 
1 5 °C, and can be treated as a function of temperature for larger temperature changes), 
10 AT is the difference between the subsequently measured temperature of the 

sample and the initial temperature, and 

dn/dC is partial derivative of index of refraction, n, of the multi-component 
liquid mixture, with respect to concentration of the given component, evaluated from 
known data in a region near the initial concentration and near the sample temperature. 
15 Also provided is a method for determining the temperature of a given 

sample of a liquid, which can be either a pure liquid or a multi-component liquid 
having a substantially constant concentration of the components. The method 
includes the step of determining the refractive index of the sample of the liquid as 
initially described above, and then comparing the refractive index of the sample of the 
20 liquid to predetermined data relating different temperatures of the liquid to 
corresponding values of the refractive index of the liquid. In this way, the 
temperature of the sample of the liquid can be determined which corresponds to the 
refractive index of the sample of the liquid which was determined as described above. 

The present invention yet further provides an apparatus for determining 
25 the concentration of a solute and the temperature of a light transmitting liquid. The 
apparatus includes a vessel which contains the liquid, a light source, and a light beam 
position sensor. The vessel has a planar entrance side and a planar exit side with a 
known angular relationship therebetween, and is immersed in known surroundings. 
The light source is located so as to cause a beam of light to impinge on the entrance 
30 side at an angle Q t with respect to the normal to the outer surface of the entrance side 
where the beam impinges, to pass through the material, and to then pass through the 
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exit side, from which it exits at an angle 8 e with respect to the normal to the outer 
surface of the exit side where the beam exits. The sensor determines the exit angle 8 e , 
thus permitting the refractive index of the material to be determined, based on the 
angles and 0 e , by applying SnelPs law. The application of SnelFs law occurs at the 

5 interface between the surroundings and the entrance side, the interface between the 
entrance side and the material, the interface between the material and the exit side, 
and the interface between the exit side and the surroundings. The apparatus includes a 
computer which is suitably programmed to calculate the refractive index of the liquid, 
and to determine the unknown concentration or temperature of the liquid sample 

1 0 based on comparison with predetermined data. 

The present invention thus provides methods and an apparatus capable 
of carrying out precision refractive index, concentration and temperature measurement 
of liquids. The apparatus and methods overcome the disadvantages of prior systems 
and methods. In particular, the present method and apparatus permit real time, non- 

1 5 invasive, remote measurements with high resolution and can be carried out using 
simple and inexpensive components. Further, the apparatus is compact and easy to 
operate. Yet further, measurements can be based on readily quantifiable and 
measurable parameters, rather than visual, subjective parameters. Real-time remote 
monitoring of liquid component concentration, such as real-time mixture monitoring, 

20 precision solute and contaminant analysis, and measurement of diffusion coefficients 
of components of a multi-component liquid mixture, for example, are all possible with 
the present invention. 

These and other features and advantages of the present invention will 
be pointed out in the following specification, taken in connection with the 

25 accompanying drawings, and the scope of the invention will be set forth in the 
appended claims. 


BRIEF DESCRIPTION OF THE DRAWINGS 
Figure 1 A is a schematic diagram of an exemplary embodiment of an 
apparatus according to the present invention; 
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Figure IB is a detailed view of a vessel of the apparatus of Figure 1 A, 
showing various parameters of interest; 

Figure 2 is a plot of refractive index of an ethanol-water mixture versus 
ethanol concentration, determined using the present invention and compared to known 
5 values; 

Figure 3 is a figure similar to Figure 2 for a methanol-water mixture; 

Figure 4 is a graph of measured concentration vs. true concentration 
for aNaCl-water solution; 

Figure 5 is a figure similar to Figure 4 for a MgCl 2 -water solution; and 
1 0 Figure 6 is a plot similar to that of Figure 5, but for very small 

concentration changes. 


DETAILED DESCRIPTION OF THE INVENTION 
Reference should now be had to Figures 1 A and IB for a description of 
the exemplary embodiments of the apparatus and method of the present invention. An 
1 5 exemplary apparatus for determining the refractive index of a material, according to 
the present invention, is designated generally by reference numeral 10. The apparatus 
includes a vessel 12 which contains a multi-component liquid 14 (e.g., a solution). As 
best seen in Figure IB, the vessel 12, which can be, for example, a precision cuvette, 
has a planar entrance side 16 and a planar exit side 1 8. The vessel 12 can be 
20 immersed in known surroundings 20, such as ordinary atmospheric air. 

Referring again to Figure 1 A, the apparatus 10 also includes a laser 
light source 22 which is located so as to cause a collimated beam of coherent light 24 
to impinge on the entrance side 1 6 at an angle 0; with respect to a normal 26 to the 
outer surface of the entrance side 1 6 where the light beam impinges. Light source 22 
25 is located with respect to vessel 12 such that beam 24 then passes through the liquid 
1 4 and then through the exit side 1 8, from which it exits at an angle 6 C with respect to 
a normal 28 to the outer surface of the exit side 1 8 where the beam exits. 

Apparatus 10 further includes a beam position sensor, designated 
generally as 30, which determines the exit angle 0 C , thus permitting the refractive 
30 index of the material 14 to be determined, based on the angles Q { and 0 C , by using 
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Snell's law. In particular, Snell's law is applied at an interface 32 between the 
surroundings 20 and the outer surface of the entrance side 16, an interface 34 between 
the inner surface of the entrance side 16 and the liquid 14, an interface 36 between the 
liquid 14 and the inner surface of the exit side 1 8, and an interface 38 between the 
5 outer surface of the exit side 1 8 and the surroundings 20. 

It should be noted that the vessel 12 need not necessarily be completely 
filled with the liquid 14, it is only necessary that liquid 14 be present in sufficient 
quantity to contain the path of the beam of light 24. Further, it will be appreciated 
that the light source can be a laser, so as to provide a collimated beam of coherent 

1 0 light 24, and that the laser may be a helium-neon laser having a wavelength of 632.8 
nm. Of course, other types of light sources emitting collimated light beams at other 
wavelengths can be employed. 

As noted, the planar entrance side 1 6 and the planar exit side 1 8 of the 
vessel 12 have a known angular relationship between them, such that Snell's law can 

15 be applied. The entrance and exit sides 16, 18 may be perpendicular, so as to make 
various calculations simpler, as will be apparent below. Sensor 30 preferably includes 
a light beam position sensor 40 which senses a position of the beam of light 24 after 
the beam exits the exit side 1 8 of the vessel 12. The beam is designated as 42 in this 
region. Position sensor 40 preferably has a known geometrical positional relationship 

20 to the vessel 12 so that the sensed position of the beam 42 can be related to the exit 
angle 8 C . In a simple form of the invention, the beam position sensor 40 can simply 
be a suitable flat surface wherein the position of illumination by the beam 42 can be 
detected and related to an initial position. In a preferred form, a semiconductor beam 
position sensor, to be discussed below, is employed. 

25 It is envisioned that the present invention will have particular utility for 

measurements of temperature and concentration to be carried out with materials 14 
which are in the liquid phase. Measurement of index of refraction of gases may be 
accomplished by containing the gases in a sealed vessel at a desired pressure (in 
practice, the near-unity index of refraction of most gases may limit the utility of the 

30 method except, possibly, at high pressures). 
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The entrance and exit sides 16, 18 of the vessel 12 preferably have a 
substantially uniform thickness, t, as shown in Figure IB, and are preferably made 
from the same light transmitting material such that they have substantially the same 
refractive index. Although the dimension t is shown on a side of the vessel other than 
5 the entrance and exit sides 1 6, 18 in Figure IB, it will appreciated that only these 
particular sides need have the uniform thickness and refractive index to obtain the 
desired benefits. However, in a precision optical cuvette, all sides are uniformly thick 
and have a substantially uniform refractive index, and it is such a cuvette which is 
envisioned as being one preferred form of vessel for use with the present invention. 
1 0 The apparatus 1 0 can also include a computer 44 which receives a 

representation of the position of the beam 42 of coherent light and which determines 
the refractive index of the liquid 14, according to the following equation: 

d(n,) = [W 0 -W( ni )] cos 8^ + {D-[W 0 -W( n| )]} tan^n,^] (2) 

where: 

1 5 d(ni) is the position of said beam of light 42, measured with the position sensor 

40, relative to an initial calibration value (which can be measured with a test liquid, as 
discussed below), 

W 0 is a reference value of distance between an inner corner 60 formed by the 
intersection of the planar entrance and exit sides 1 6, 1 8, and a point where the beam 
20 of light 42 f exits the exit side 1 8 (see interface 38'), determined by measurement and 
calculation, as discussed below, 

6 e o is a reference value of the angle at which the beam of light 42' exits the 
exit side 18, determined during calibration, 

D is distance from the point where the beam of light exits the exit side to the 
25 position sensor 40, measured during calibration, 

W(n,) is a value of the distance between the inner corner 60 formed by the 
intersection of the entrance and exit sides 16, 18, and a point where the beam of 
coherent light 42 exits the exit side 1 8 when the vessel 12 contains the liquid being 
measured, during a given measurement, defined according to the equation: 
30 W(n ( ) = L cot[arc sin((sin B^nJ] + 1 tan{arc sin[(n/ig(l - (sin O^n,)*]} (3) 

where: 
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L = distance from the inner comer 60 to a point where the beam of light passes 
from the entrance side into the liquid (interface 34), determined by measurement and 
Snell's-law calculation, 

n, = refractive index of liquid which is to be determined, 
5 t = thickness of vessel walls, and 

n g ~ refractive index of vessel walls, 

and 

6 e (n,) is the value of 0 e , during measurement of a given liquid, defined 
according to the equation: 

10 Oefa) = arc sin {n, cos[arc sin ((sin 8/^]}. (3) 

It will be appreciated that equations (2), (3) and (4) hold for both 
positive and negative values of d(n,). Further, it will be appreciated that a neat, 
explicit solution for d(n,) is not available so that equations (2), (3) and (4) can be 
solved, for example, by manual or computer-aided iterative calculation. 

1 5 The position sensor 40 is preferably a semiconductor light beam 

position sensor, such as a photodiode beam position sensor to be further described 
below, coupled to a position-to-voltage converter 46 which produces an output 
voltage which is substantially linearly proportional to the displacement of the beam 42 
from its reference position 42' determined during calibration. Further, computer 44 is 

20 preferably a general purpose digital computer, although those of skill in the art will 
appreciate that an application specific integrated circuit (ASIC) can also be employed. 
Apparatus 10 can further include a first digital voltmeter 48 which receives the output 
voltage from the position-to-voltage converter 46 and which supplies a corresponding 
digital signal to the computer 44 as the aforementioned representation of the position 

25 of the beam of light 42. 

Apparatus 10 can further include a temperature sensor 50 which is 
positioned to sense the temperature of the liquid 14 and which is coupled to the 
computer 44 so as to supply the computer 44 with a digital signal representative of the 
temperature of the liquid 14. The coupling between the temperature sensor 50 and the 
30 computer 44 can be, for example, a second digital voltmeter 52. Sensor 50 is 
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preferably a thermistor, but could also be a thermocouple, or any other type of 
suitable temperature sensor known to those of skill in the art. 

Apparatus 10 can also include a rotation stage 54 on which the vessel 
12 is mounted. Rotation stage 54 permits the angle of incidence of the beam of light 

5 42 onto the vessel 12 to be adjusted during calibration of the apparatus 1 0, such that 
the reference position of the beam from the exit side 18 of the vessel 12 during 
calibration is substantially perpendicular to the plane of the beam position sensor 40 
and impinges the sensor 40 at its midpoint. The beam 42' shown in Figures 1 A and 
IB represents a reference beam during a calibration process, which has been adjusted 

10 to be substantially perpendicular to the beam position sensor 40, and will be described 
more fully below. 

For purposes of validation and calibration, apparatus 10 can further 
include a scale 58 interposed between the rotation stage 54 and a mounting surface 56. 
Scale 58 is preferably a digital balance. It can be provided in order to determine the 

15 mass of the liquid 14 being measured, and particularly changes in mass concentration 
of a liquid being measured, by weighing samples of the liquid 14 contained in the 
vessel 12 as additional mass of a given component is added to change the 
concentration. Again, the scale is not normally envisioned for use in a production 
environment, since the apparatus of the present invention would be employed to 

20 determine concentration. 

It should be noted that concentration of a component in a multi- 
component liquid may be expressed in molar, volume percent or mass percent units in 
accordance with the present invention, however, mass percent is considered to be 
preferable. 

25 The apparatus of the present invention can be employed to determine 

the unknown temperature of a sample of liquid 14 by measuring its index of 
refraction. The index of refraction of the material 14 can then be compared to 
predetermined data which relate different temperatures of the material to 
corresponding values of the refractive index of the material. Thus, the temperature, 

30 hitherto unknown, can be determined for a given sample of liquid by measuring the 
refractive index of the sample and converting the measured refractive index to a 
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corresponding sample temperature. Computer 44 can be programmed to carry out 
such calculations, although they can also be done manually. Temperature 
determination is most practical for pure materials, or for multi-component materials 
having a substantially constant concentration of the components. The refractive index 
5 of the unknown material is measured, and is then compared to known values of 
temperature vs. index of refraction, for example, using interpolation, a curve fit, or 
any other suitable means which will be apparent to those of skill in the art. Straight 
line, polynomial or spline curve fits can be employed. The "Tablecurve 2D" program 
available from SPSS software has been found to be useful for such curve fitting. 

1 0 The apparatus 1 0 according to the present invention can also be used to 

determine the unknown concentration of a given component of a multi-component 
liquid mixture. In this case, the liquid 14 is a multi-component liquid mixture and the 
determination is carried out by comparing the measured refractive index of the multi- 
component liquid mixture to predetermined data which relate different concentrations 

15 of the given component of the multi-component liquid mixture to corresponding 

values of the refractive index. Again, this comparison can be carried out manually, or 
preferably by a suitably programmed computer 44. Thus, an unknown value of the 
concentration of a given component of a sample of the multi-component liquid 
mixture can be determined by measuring the refractive index of the sample and 

20 converting the measured refractive index to a corresponding concentration of the 
given component in the sample of the multi-component liquid mixture. Again, as for 
temperature, a suitable interpolation between tabulated values, comparison with a 
curve fit, or any other desired method can be employed. 

As will be discussed more fully below, the refractive index of a multi- 

25 component liquid mixture is a function of the mixture temperature, mixture 

component concentration, and wavelength of the incident light. While the wavelength 
of the incident light is substantially constant when using a laser, the temperature and 
concentration can vary. When it is desired to determine the concentration of a given 
component in the mixture, in the presence of temperature changes, suitable correction 

30 must be carried out For relatively large changes in concentration, values of refractive 
index as a function of concentration can be stored in computer 44 (or otherwise) for a 
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number of different values of temperature. The measured refractive index of the 
specimen of unknown concentration can then be compared to predetermined data for 
at least two temperatures nearest the temperature of the specimen, which can be 
measured with the temperature sensor 50, Interpolation can then be carried out 
5 between the values of the predetermined data at the at least two nearest temperatures 
so as to obtain the value of the concentration of the given component of the liquid 14 
for the actual temperature of the liquid 14. An ordinary bilinear interpolation, familiar 
to thermodynamics students who have used the steam tables for the properties of 
water vapor, can be employed and can be programmed into computer 44, as will 

1 0 readily be appreciated by those of skill in the art. Other suitable interpolation 
methods employing curve fitting and the like can also be used. 

Reference should now be had to Figure 3, which will be discussed in 
greater detail below. Note that in some liquid mixtures, for example, the methanol- 
water mixture shown in Figure 3, a given value of the refractive index of the mixture 

1 5 can correspond to two different values of concentration. In this case, if no other 

information is known, a unique value of concentration for a given a measured value of 
refractive index cannot be determined. However, if the previously measured values of 
concentration, which indicates a range of expected concentrations are available, then 
the correct value of the concentration can be selected. For example, computer 44 

20 could be programmed to store values of the concentration of the given component of 
the liquid taken at previous proximate points in time and to employ the previous 
values to select the correct value of the concentration when the concentration is 
double valued with respect to the index of refraction. 

Still referring to Figure 3, for example, consider the dark dots which 

25 represent a first series of measurements carried out from 0 mass % methanol through 
60 mass % methanol. Suppose, for an industrial process, it was known that the initial 
refractive index was about 1 .33 1 for an initial methanol concentration of 0 mass %. 
Further, suppose that measurements of indices of refraction were made every 10 
seconds, and the next two measurements yielded refractive indices of approximately 

30 1 .332 and 1 .333, which according to the data in Figure 3 could correspond to about 5 
and 10 mass % methanol, respectively, or to about 84 and 83 mass % methanol, 
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respectively. Assuming that it is known that the concentration of methanol increases 
relatively slowly with time (and in the absence of radical jumps or piecewise 
discontinuities in the relationship between n and C), it may be inferred that 5 and 10 
mass % are the correct methanol concentrations corresponding to the refractive index 

5 measurements of 1 .332 and 1 .333, respectively. When a value of 1 .335 was measured 
as the refractive index approximately 10 seconds later, based on the previous data, it 
may be inferred that this corresponded to a 1 5 mass % concentration of methanol 
(instead of 82 mass %), according to the data in Figure 3. Other methods of dealing 
with double-valued systems are possible, such as deliberately introducing a small 

1 0 increase in concentration of a given component (e.g., methanol) and measuring a 

change in the index of refraction. An increase means the left side of curve in Figure 3 
while decrease means the right side of curve. A decrease in concentration could also 
be introduced to determine which side of the n vs. C curve should be used to indicate 
the concentration corresponding to a measured value of refractive index. Thus, the 

1 5 slope of the n vs. C curve is used to resolve the double - valued problem. 

As noted, the apparatus 10 can be used for measurement of a multi- 
component liquid 14. Further, apparatus 10 can also include the aforementioned 
computer 44 and the aforementioned temperature sensor 50 which, as noted, can be 
positioned to sense temperature of the multi-component liquid 14 and which can be 

20 coupled to the computer 44 so as to supply the computer with a signal which is 
representative of the temperature of the multi-component liquid 14. When it is 
expected that relatively small concentration changes will be encountered, the 
apparatus 10 can be used to calculate changes in the concentration of given 
components of the multi-component liquid mixture from an initial value. Computer 

25 44 can be programmed to determine the refractive index of a sample of the multi- 
component liquid mixture based on the entrance and exit angles 0, and 0 e , by applying 
Snell's law at the interfaces 32, 34, 36, 38 as discussed above. Computer 44 can 
further be programmed to determine a change in the refractive index of the multi- 
component liquid mixture by subtracting, from the refractive index determined by the 

30 computer, an initial refractive index of the sample of the multi-component liquid 
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mixture corresponding to an initial concentration of a given component at a given 
temperature measured with the temperature sensor 50. 

Computer 44 can be further programmed to receive the signal which 
represents the temperature of the given sample of the multi-component liquid mixture, 

5 and to determine the change, AC, in the concentration of the given component of the 
sample of the multi-component liquid mixture from the initial concentration according 
to Equation (1) above. Equation (1) is advantageously used in cases where AC is 
relatively small. As a rule of thumb, Equation (1) is appropriately used for up to 5 or 
even 10% changes in C, although for substantially constant values of dn/dc (at a 

10 constant temperature), Equation (1) can be used for almost any value of AC. 

In view of the foregoing description, it will be appreciated that the 
apparatus 1 0 can be employed to measure the index of refraction of a liquid, the 
temperature of a liquid, the concentration of a given component in a multi-component 
liquid, or a change in concentration of a given component in a multi-component 

15 liquid. In all cases where various calculations are referred to as being performed by 
computer 44, it will be appreciated that they could also be performed manually or by 
any other known means. 

A method, according to the present invention, for determining the 
refractive index of a liquid 14, will now be discussed. The method includes the steps 

20 of (a) causing the material to be contained in a vessel 12 having an entrance side 1 6 
and an exit side 1 8 of known refractive index; (b) causing a beam of coherent light 24 
to impinge on the entrance side 1 6 at an angle 0j with respect to a normal 26 to the 
surface of the entrance side 16 where the beam impinges; and (c) calculating the 
refractive index of the liquid 14. 

25 In step (a), as for the apparatus 10, the vessel 12 can have entrance and 

exit sides 16, 1 8 with a known angular relationship therebetween and can be 
immersed in known surroundings 20. In step (b), the beam of light 24 can pass 
through the liquid 14, and then through the exit side 1 8 of the vessel 12, from which it 
exits at an angle 0 e with respect to a normal 28 to the surface of the exit side 1 8 where 

30 the beam exits. The angle 0 C can be measured directly or indirectly within the scope 
of the present invention. The refractive index of the material can then be calculated 
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based on the angles and 6 C , by applying SnelFs law at the interfaces 32, 34, 36, 38, 
as described discussed above. 

In a case where the liquid 14 is a multi-component liquid mixture, a 
further step (d) can be performed, which comprises repeating steps (a) through (c) for 
5 samples of the multi-component liquid mixture having different concentrations of a 
given component. Accordingly, a curve of values of the refractive index of the multi- 
component liquid mixture corresponding to different concentrations of a given 
component can be obtained. Such curves are shown, for example, in Figures 2 and 3 
and will be discussed further below. Throughout this application, "curve" should be 

1 0 understood broadly to include a physical, graphical plot, tabulated data values, 
parameters stored in a computer, and the like. Yet further, step (d) can be repeated 
under substantially isothermal conditions at a plurality of different temperatures, so as 
to obtain a family of curves of values of the refractive index of the multi-component 
liquid mixture corresponding to different concentrations of a given component for the 

1 5 plurality of different temperatures. For example, a family of curves similar to Figure 
2 could be obtained at temperatures other than 20°C. The data obtained by this 
method could then be used, for example, by computer 44, to determine the unknown 
concentration of a given sample by measuring its index of refraction and, optionally, 
temperature. 

20 In another aspect of the invention, the liquid 14 can be either a pure 

liquid or a multi-component liquid mixture having a substantially constant 
concentration of its components. An additional step can be performed in this case, of 
repeating steps (a) through (c) for samples of the liquid at different temperatures, so as 
to obtain a curve of values of the refractive index of the liquid corresponding to 

25 different temperatures of the liquid. While such a curve is not shown in the drawing 
figures herein, it will represent a plot of refractive index versus temperature for either 
a pure liquid or a multi-component liquid mixture where the concentrations of all the 
components are constant. Further discussion on the refractive index as a function of 
temperature, concentration and wavelength will be presented below. However, as 

30 used herein, "substantially constant concentration" means that any changes in 
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concentration of the components are small enough such that they may be neglected in 
determining a change in the refractive index. 

The method can also include the step of providing an beam position 
sensor 40 as discussed above, which senses a position of the beam of light 42 
following exit from the exit side 1 8. Just as for the vessel 12, the entrance and exit 
sides 16, 18 can be substantially perpendicular. Sensed position from the beam 
position sensor 40 can be used in calculating the refractive index. 

As noted with respect to the vessel 12, the entrance and exit sides 1 6, 
18 of the vessel 12 can have a substantially uniform thickness, such as t, and a 
substantially uniform refractive index. Step (c) can include detennining the refractive 
index of the liquid sample 14 according to equations (2), (3) and (4) set forth above. 

It is desirable that apparatus 1 0 be calibrated before use, and one 
possible method of performing such calibration will now be discussed. Values of the 
parameters W 0 and 9^ can be determined. The calibration can include the steps of 
causing a first liquid of known refractive index to be contained in the vessel 12. A 
relative rotation between the vessel 12 and the position sensor 40 can then be carried 
out, for example, using the rotation stage 54, such that, for the first liquid, the beam of 
light 42 exits the exit side 18 substantially perpendicular to the beam position sensor 
40 and impinges thereon at its midpoint Accordingly, the aforementioned initial 
calibration value of d may be taken as 0. 

The calibration method can further include recording W 0 and 0 eO , 
emptying the vessel 12 of the first liquid, and causing a second liquid of known 
refractive index to be contained in the vessel 12. The refractive index of the second 
liquid is different than that of the first liquid. Values of d, W and 8 e corresponding to 
the second liquid can then be determined, with the values of W and 9 e being 
respectively calculated from the equations (3) and (4) above. The temperatures of the 
first and second liquids can be measured, during calibration, using temperature sensor 
30 such that the refractive index of each liquid is accurately determined at well- 
defined temperatures. Calibration can be carried out under constant temperature 
conditions. 
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The parameters L and W 0 can be determined as follows. The value of 
L can be determined by measuring the position L'of the laser beam 24 on the outside 
of the vessel 12 and then using Snell's law in conjunction with the vessel wall 
thickness t to determine the value of L as shown in Figure IB. One subtracts from L' 
the downward displacement of the beam through the entrance side 16, determined 
from Snell's law. Once L is known, Snell's law can again be used to determine W 0 . 
Alternatively, the wall thickness t can be subtracted from W 0 ' to yield W„. 
Accordingly, the values of L and W 0 may be determined by measurements and 
calculations using Snell's law. The measurements of W 0 ' and L' can be made with a 
ruler. W 0 ' and L 1 can be measured from an actual corner of a "virtual" sharp corner 
where the planar sides would intersect when rounded corners are present. 

The present invention also provides a method for determining the 
concentration of a given component in a sample of a multi-component liquid mixture. 
The method includes the steps of determining the refractive index of the sample of the 
multi-component liquid mixture as described above, and then comparing the refractive 
index of the sample thus determined to pnrietermined data relating different 
concentrations of the given component of the multi-component liquid mixture to 
corresponding values of the refractive index of the multi-component mixture, so as to 
convert the determined refractive index of the sample to a corresponding 
concentration of the given component in the sample of the multi-component liquid 
mixture. It will be appreciated that the above-described apparatus of figures 1A and 
IB can be used in carrying out this method. The aforementioned comments on 
calibration also apply to this method. 

This method can also include the additional step of recording the 
temperature of the sample of the multi-component liquid mixture. In this latter case, 
the step of comparing the refractive index of the sample of the multi-component 
liquid mixture to the predetermined data can include comparing the refractive index of 
the sample previously determined to values of the predetermined data relating the 
different concentrations of the given component to the corresponding values of the 
refractive index for at least two temperatures near the recorded temperature of the 
sample, and then interpolating between those values. As described above with respect 
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to the apparatus 10, the interpolation can be carried out at the at least two 
temperatures near the recorded temperature of the sample of the multi-component 
liquid mixture, so as to obtain the value of concentration of the given component in 
the sample of the multi-component liquid mixture which is determined for the 
5 recorded temperature. Again, as described above, the temperature can be determined 
with the temperature sensor 50. 

In the foregoing method, the predetermined data relating different 
concentrations of a given component of the multi-component liquid mixture to 
corresponding values of the refractive index of the multi-component liquid mixture 

1 0 may be double valued as described above, such that certain values of the refractive 
index correspond to two different values of concentration. For example, with 
reference to Figure 3, refractive indices from about 1.331 to about 1.342 each 
correspond to two different values of methanol concentration. In this case, the 
method can further include storing previously measured values of concentration for 

1 5 the given component in the sample of the multi-component liquid mixture taken at 
proximate prior points in time, and then using such stored previously measured values 
of the concentration to indicate an expected range of a newly measured concentration 
value so that the correct newly measured concentration value can be selected. 

The present invention further provides a method for determining a 

20 change in concentration of a given component in a sample of a multi-component 
liquid mixture from an initial concentration corresponding to an initial refractive 
index of the sample at an initial temperature. This method is advantageous for 
relatively small changes in concentration, as discussed above with respect to Equation 
(1). This method includes the steps of (a) determining the refractive index of the 

25 sample of the multi-component liquid mixture as described above; (b) determining a 
change in the refractive index, An, by subtracting the refractive index determined in 
the preceding step from an initial refractive index; (c) measuring the temperature of 
the sample of the multi-component liquid mixture, for example, with the temperature 
sensor 50; and (d) determining the change, AC, in the concentration of the given 

30 component in the sample of the multi-component liquid mixture, from the initial 
concentration, according to equation (1) given above. 
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The present invention yet further provides a method for determining 
the temperature of a sample of a liquid, wherein the liquid may be a pure liquid or a 
multi-component liquid having a substantially constant concentration of each 
component. The method includes the steps of (a) determining the refractive index of 
5 the sample of the liquid as described above; and comparing the determined refractive 
index of the sample to predetermined data which relate the different temperatures of 
the liquid to corresponding values of the refractive index of the liquid, so as to convert 
the determined refractive index to a corresponding temperature of the sample of the 
liquid. Again, the above-described apparatus 10 can be employed in carrying out this 
10 method. The above-described equations (2), (3) and (4) can be employed. The 
foregoing description of calibration of the apparatus 10 is also applicable to this 
method. 

Further description of the theory of the present invention will now be 
set forth. Both absolute and relative index of refraction changes can be measured with 

15 the present invention. Extremely small changes in concentration can be resolved, on 
the order of 0.01 mass %, while at the same time precise measurements can be made 
over a large range of concentration values (0-100 mass %). The apparatus of the 
present invention can be readily adapted for continuous flow measurement, making it 
suitable for real-time monitoring and control of liquid manufacturing processes. 

20 Concentrations may be mass percent, volume percent, or molar; mass 

percent is preferred. Mass concentration may be known, e.g., by measuring the mass 
of liquid(s) added to the vessel 12 with scale 58, so that the refractive index as a 
function of concentration is obtained. Conversely, if the relationship between 
concentration and refractive index is already known, the concentration of the mixture 

25 can be determined. The refractive index of a liquid depends on its density, as 
described in J. D. Spear, R. E. Russo, and R. J. Silva, "Collinear Photothermal 
Deflection Spectroscopy with Light-Scattering Samples," 29 Applied Optics 4225-34 
(1990), and the wavelength of the incident light, as set forth in G. M. Hale and M. R. 
Querry, "Optical Constants of Water in the 200-nm to 200-|im Wavelength Region," 

30 12 A pplied Optics 555-63 (1973). Fluctuations in both temperature and concentration 
will change the liquid density (with pressure effects neglected). Thus, for an 
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isothermal, achromatic measurement, the change of the refractive index of a multi- 
component liquid is determined only by the concentration of the components. 

In general, the measured refractive index of a multi-component liquid 
mixture is a function of its temperature, T, concentration of a given component, C, 
5 and the wavelength of the incident light, X: 

n = n(T,C,A) (5) 
From Eq. 5, the change of refractive index, Aw, of a multi-component liquid mixture 
is: 

An * (fln/aT) AT + (dn/3C)AC + (dn/dX)AX (6) 

10 for small changes in the temperature, concentration, and wavelength. Since laser light 
is highly monochromatic and wavelength-stable, AA~0, and wavelength effects can be 
neglected. As noted, the concentration herein is preferably defined on a mass percent 
basis. Thus, for liquid components A and B: 

C A = m A /(m A + m B ) and C B = m B /(m A +m B ), (7ajb) 

1 5 where C A , C B are concentration of components A and B, respectively, and m A and m B 
are masses of components A and B, respectively. Thus, by measuring the change in 
refractive index and the liquid temperature, the concentration change can be expressed 
using Eq. (4) above. 

To measure Ait, SnelFs law is employed to relate the incident and exit 

20 light beam angles, and the liquid refractive index as the beam passes through the 
vessel 12. For light passing from medium 1 to 2 with indices of refraction n } and « 2 , 
respectively, SnelFs law takes the form 

sin 6 X = «2 sin 0 2 (8) 
where 6t and 6 2 are the angles of the light beam in mediums 1 and 2, respectively, 

25 measured with respect to the normal to the interface of medium 1 and medium 2. 
In Figure IB, the laser beam 24 is sent into the vessel 12 at an incident angle 0 { such 
that it passes through the wall 16 into the material 14 and then strikes and passes 
through the wall 1 8, which may be perpendicular to the entrance wall 1 6. Applying 
SnelFs law to this optical path, the exit position Wand exit angle 0 C are given by 

30 equations (3) and (4) above. 


WO 00/26665 PCT/US98/23179 

22 

Note that both 0 e and W are (in general) required, because a change in 
n, results in a change both in the exit angle 8 C from the vessel 12 and a translation of 
the beam 42 along the exit wall of the vessel 12. Both will result in a beam position 
change at the beam position sensor 40. The change in the beam position, d, at the 
5 sensor 40 in Figure la can be expressed by equation (2) above. 

In many cases, 6 e is more significant than W, and a simplified, explicit 
formula for the unknown refractive index n, can be used: 

n, = (sin 2 8i + sin 2 8 e ) w . (9) 
The importance of lateral beam displacement on the vessel exit side 1 8 
1 0 (difference between W and W 0 ) decreases with increasing distance D between the 
vessel 12 and beam position sensor 40. For example, with D = 2 cm, lateral beam 
displacement contributes about 4% of total beam deflection; this declines to about 1% 
at D - 10 cm and continues to decline with increasing D. 

Referring to Eq. (1), temperature changes in the liquid will result in a 
1 5 perceived concentration change. Rather than attempt to maintain an isothermal 

system, the liquid temperature can be monitored, and the concentration change can be 
compensated for temperature by using Eq. (1) and values of dn/dT such as those in 
Table 1 below. It will be seen that values of total derivative, dn/dT, are tabulated in 
Table 1 because the data therein is for pure liquids. Values for mixtures can be 
20 estimated by linearly interpolating, using the values of dn/dT for each component. 
The value of dn/3T may often be approximately constant with concentration, or can be 
treated as a function of concentration and a double interpolation carried out 


Table 1. Thermo-Optical Properties of Liquids 



n(589nm,20°C) 

n(632.8 nm, 20 °C) 

10 4 x dn/dT (K 1 ) 

water 

1.3330 [17] 

1.331 [21] 

-1.04 [18] 

ethanol 

1.3614 [17] 

1.358 [21] 

-3.9 [21] 

methanol 

1.3290 [17] 

1.325 [21] 

-3.9 [21] 


EXAMPLE 1 

A Melles-Griot 05-LHR-l 1 1 1 mW linearly polarized HeNe laser with 
a wavelength of 632.8 nm served as the light source 22. A precision square glass 
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cuvette, with interior dimensions of 10 x 10 x 35 mm, was used as the vessel 12 (and 
optical cell). The lengths of each cuvette side were measured and found to be 
nominally 10mm with an uncertainty of ±25 urn, and the edges were square to within 
±3.0 mrad. The laser beam 24 was directed onto the cuvette wall 1 6 at a specified 
5 incident angle ranging from 70 to 80° (measure with respect to the normal to the wall) 
using a rotation stage 54. The beam 42 exits the adjacent cuvette wall 1 8 at an angle 
0 e that depends on the index of refraction, n„ of the liquid in the cuvette, and the 
exiting beam strikes a semiconductor beam position sensor 40. 

As the concentration - and hence the refractive index - of the liquid in 

1 0 the cuvette changes, the beam position on the sensor 40 will move accordingly. A 
UDT LSC-5L linear beam position sensor with a 1 5 V bias voltage was used as sensor 
40 to monitor this position change of the beam 42. This highly linear, low noise 
position sensor was of the single photodiode type and was capable of detecting beam 
displacements as small as 2 urn. It was manufactured by United Detector 

1 5 Technologies, Inc. of San Jose, California. The sensor output was sent to a position- 
to-voltage converter 46, which produced an output voltage linearly proportional to the 
beam position. This voltage output was in turn measured using a 6.5 digit Keithley 
Model 2000 digital voltmeter (DVM), shown in Figure 1 A as first voltmeter 48. The 
measured sensor sensitivity was 3.89 mV/um. The DVM readings were sent via a 

20 Hewlett-Packard industry standard general purpose bus (GPIB) interface to computer 
44, which was a personal computer running Visual Basic for data acquisition. 
Typically, 10 readings were taken over a 1 sec interval and averaged to obtain a final 
reading. 

The UDT SL5-2 beam position sensor has three leads. The center lead 
25 can be grounded while the other to leads A and B each provide a current I A and I B , 
respectively. With the beam 42 centered on the beam position sensor 40, the same 
current (I A = I 8 ) will flow in each of the leads A and B. A change in the currents will 
occur when the beam 42 goes off center. The currents I A and I B can be converted to 
respective voltages V A and V B , respectively, using a well known operational amplifier 
30 technique. The differential voltage, V A - V B , directly relates to the beam position. To 
avoid errors in measuring caused by fluctuations in light beam intensity, the voltage 
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difference (V A - V B ) is normalized by the voltage sum (V A + V B ). This can be 
accomplished by adding V A and V B by using a well known operational amplifier 
technique and by using an AD734 4-quadrant multiplier-divider manufactured by 
Analog Devices, Inc. of Norwood, Mass. to provide a quantity (V A -V B )(V A +V B) , 
5 which is also directly proportional to beam position relative to the center of the LSC- 
5L beam position sensor. 

The concentration of the given component in the liquid (mass % 
concentrations were used) was determined using an Ohaus laboratory digital balance, 
shown in Figure 1A as scale 54, with a resolution of 0.01 g. The glass cuvette 12 and 

1 0 rotation stage 54 were mounted on scale 58. Very small concentration changes that 
would otherwise not be resolvable with the scale 58 were made by first preparing a 
diluted mixture of known concentration, and then adding this diluted mixture to the 
cuvette, so that the total mass was large enough to be measured by the scale 58. To 
minimize vibration effects, the system was assembled on a vibration-isolated optical 

1 5 bench. Vibration isolation is believed to be desirable in all embodiments of the 
present invention. A 2 mm diameter thermistor was used as temperature sensor 50 
and was placed under the test liquid surface inside the cuvette to monitor the liquid 
temperature. The liquid temperature was recorded simultaneously with the beam 
position using a second Keithley multimeter, shown in Figure 1 A as voltmeter 52. 

20 To calibrate the system, two test liquids, pure water and methanol, 

which both have well-characterized refractive indices, were used. Calibration needs 
to be done only once, provided the components in the optical path of the system are 
not moved. The pure water was first placed in the vessel 12, and the beam position 
sensor 40 was oriented perpendicular to the exit laser beam 42, which impinged on the 

25 center of the sensor 40. The beam position W 0 and exit angle 6^ required in Eq. (2) 
were also recorded. Next D in Eq. (2) was determined by removing the water and 
adding pure methanol to the cell, which resulted in a beam displacement, d, at the 
beam position sensor 40. The new values of W and 0 e were determined from Eqs. (3) 
and (4), which allowed for determining the sole unknown, D. The liquid temperature 

30 during calibration was also measured to ensure accurate values of n. 
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To demonstrate the refractive index measurement capabilities of the 
present invention, ethanol-water and methanol-water mixtures were used, as they have 
refractive indices that exhibit a strong non-linear dependence on the concentration, as 
described in D. J. Williams, Nonlinear Optical Properties of Organic and Polymeric 

5 Materials (American Chemical Society, Washington, D.C., 1983), and as is well 
documented in the literature. For the concentration measurement (described in 
Example 2 below), two aqueous salt solutions, NaCl-H 2 0 and MgCl 2 -H 2 0, were used, 
because dn/dC for such mixtures is substantially constant, which simplifies the 
analysis and allows for determining the smallest possible resolution of the system. 

1 0 Results of measured n, for ethanol-water and methanol-water are 

shown in Figure 2 and Figure 3, respectively. Typical liquid temperatures ranged 
from 21-23 °C, and the solid-line data values have been temperature compensated to 
20 °C. The concentration ranges from 0 to 100% alcohol of each liquid. Two separate 
sets of data from two runs appear: the first ranges from 0 to 60% alcohol (solid dots), 

15 the second from 40% to 100% alcohol (open triangles). Because concentration 
changes were effected by adding one liquid to another, two runs were required to 
obtain data that span the entire 0-1 00% concentration range. In one run, alcohol was 
added to water, starting with pure water, while in the other run, water was added to 
alcohol, starting with pure alcohol. It is noted that in both Figure 2 and Figure 3 there 

20 is very good agreement in the data from the two runs in the overlap region between 40 
and 60%. The solid line is a measurement from the literature for the same liquid 
system, but at a wavelength of 589 nm, as described in D. R. Lide, Ed., CRC 
Handbook of Chemistry and Physics (CRC Press, Boca Raton, 1998), where the index 
of refraction is slightly larger. Note the agreement in the shape of the measured data 

25 with that of the published data in both figures is good. The entire concentration range 
(0-100%) was measured without any configuration changes or other adjustments to 
the system. 

EXAMPLE 2 

To measure concentration, aqueous salt solutions of NaCl and MgCl 2 
30 were used. Such solutions have a substantially linear relationship between n and C, 
i.e., dn/dC is substantially constant, as described in D. R. Lide, Ed., CRC Handbook 
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of Chemistry and Phvsics (CRC Press, Boca Raton, 1998). The values of dn/dC for 
NaCl-H 2 0 and MgCl 2 -H 2 0 solutions at 20°C are 1.71 x 10' 3 and 2.60 x 10 3 %, 
respectively. Measurements on aqueous salt mixtures permit very small concentration 
changes to be made easily, and without effects such as evaporation losses and cooling 

5 that the alcohol mixtures would incur. Known concentrations were prepared by 
measuring the salt on a microbalance and dissolving in a measured mass of water. 
Results are shown in Figure 4 and Figure 5 for NaCl-H 2 0 and MgCl 2 -H 2 0 solutions, 
respectively. To illustrate the resolution of the technique in accordance with the 
present invention, the concentration was varied by only 1%. The results are in 

1 0 excellent agreement with expected values. To determine the smallest possible 
resolution of the technique according to the present invention, tests were run on 
0.01% and 0.02% MgCl 2 -H 2 0 solutions. The results are shown in Figure 6, and it can 
be seen that the technique according to the present invention is capable of measuring 
concentration changes on the order of 0.01%. All results from the technique 

1 5 according to the present invention are in very good agreement with the true 
concentration. 

Uncertainties in the refractive index measurement arise from 
systematic errors in distance and angle measurements, and random errors such as 
vibration, air disturbances, and laser fluctuations. Random errors can be accounted 

20 for using statistical averaging techniques during data acquisition. Systematic errors 
can be determined from the following error analysis. Referring to Figures 1 A and IB, 
and Eqs. (2) - (4), uncertainties, in L, 0 i5 and D contribute to uncertainly in n,. Error 
arising from the vessel wall thickness, t, vessel index of refraction (glass was used in 
the Examples described herein), rig, and the refractive index of the two calibration 

25 liquids (water and methanol) are generally negligible. Uncertainties in L and 6j were 
taken as half of the measurement instrument resolution, and are ±0.5 mm and ±0.5 °, 
respectively. The uncertainty of the vessel-sensor distance, D, was also found to be 
±0.5 mm based on the calibration procedure discussed above. 

Uncertainty in W, U w can be obtained from Eq. (3) using the root- 

30 sum-square method: 

u w = {[ow/aL)uj 2 + [Ow/ae^Uei] 2 + [(aw/an t )u nl ] 2 >* <io) 
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Similar expressions for the uncertainty in and d can be obtained from Eqs. (4) and 
(2), respectively. The complexity of Eqs. (2) - (4) make an explicit expression for the 
uncertainty in n, difficult. Thus, the uncertainty in n-> was numerically solved for to 
derive the measured uncertainty in W and 8,. It is these results that are shown as the 
5 error bars in Figure 2 and Figure 3. The error bars in Figures 4, 5 and 6 are based on 
estimated wort-case uncertainties, determined as set forth above. Actual agreement is 
found to be better than the predicted error bars. 

When determining the concentration by measuring the index of 
refraction, the relationship between ni and C is required. The uncertainty in C, U c , is 
10 then 

u^can/acy'OJj (ii) 

Due to a larger value of dn/dC, the MgCl 2 -H 2 0 results have an uncertainty of 0.09% 
in Figure 5, which is smaller than 0.14% for the NaCl-H 2 0 solution in Figure 4. 

The thermistor used as temperature sensor 50 was calibrated in a 

1 5 separate temperature-controlled water bath using a NIST-traceable platinum- 
resistance RTD standard thermometer with a reported accuracy of 0.01 °C, resulting 
in a thermistor uncertainty of less than 0. 1 °C. The power of the HeNe laser used as 
light source 22 was 1 mW and the absorption coefficient of the test liquids, a, is on 
the order of 10* 3 , as described in D. Solimini, "Loss Measurement of Organic 

20 Materials at 6328 A," 37 Journal of Applied Physics 3314-15 (1966). Thus the 

contribution to heating by the laser beam is negligible. Temperature variations in the 
liquid were recorded and used to compensate for temperature drift of the test liquid 
during the measurement. 

Although the present invention has been described with reference to 

25 specific exemplary embodiments, it should be understood that various changes, 
substitutions and alterations can be made to the disclosed embodiments without 
departing from the spirit and scope of the invention as defined by the appended 
claims. 
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CLAIMS 

1 . A method for deteraiining the concentration of a given component in a sample 
of a multi-component liquid mixture, said method comprising the steps of: 

(a) determining the refractive index of the sample of the multi-component liquid 
5 mixture by: 

(a-1) causing the sample of the liquid mixture to be contained in a vessel 
having a planar entrance side and a planar exit side, with a known angular relationship 
therebetween, said vessel being immersed in known surroundings; 

(a-2) causing a beam of light to impinge on said entrance side at an angle 8; 
1 0 with respect to normal to an outer surface of said entrance side, to pass through the 
sample of the multi-component liquid mixture, and to then pass through said exit side, 
from which said beam exits at an angle 0 e with respect to a normal to an outer surface 
of said exit side; and 

(a-3) calculating the refractive index of the sample of the multi-component 
1 5 liquid mixture, based on said angles 8; and 0 C , by applying SnelPs law at: 

an interface between the surroundings and said entrance side; 
an interface between said entrance side and the sample of the 
multi-component liquid mixture; 

an interface between the sample of the multi-component liquid 
20 mixture and said exit side; and 

an interface between said exit side and the surroundings; 

and 

(b) comparing the refractive index of the sample of the multi-component liquid 
mixture determined in step (a) to predetermined data relating different concentrations 

25 of the given component of the multi-component liquid mixture to corresponding 
values of the refractive index of the multi-component liquid mixture, so as to 
determine the concentration of the given component in the sample of the multi- 
component liquid mixture corresponding to the refractive index of the sample 
determined in step (a). 
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2. The method of Claim 1 , wherein: 

step (a) further comprises the additional step of: 

(a-4) providing a beam position sensor which senses a position of said 
beam of light after said beam exits said exit side, and wherein step (a-3) uses 
5 said sensed position of said beam in calculating the refractive index of the 

sample of the multi-component liquid mixture. 

3. The method of Claim 2, wherein: 

said entrance and exit sides of said vessel have a substantially uniform 
thickness and a substantially uniform refractive index; and 
1 0 step (a-3) comprises determining the refractive index of the given sample of 

the multi-component liquid mixture according to a first equation: 
d(n.) = [W 0 -W(n,)] cos 8^ + {D-[W 0 -W(n,)]} tan^n,)-^], 
where: 

d(nj) is the position where said beam of light impinge on said beam position 
1 5 sensor, relative to a reference position, 

W 0 is a reference value of distance between an inner corner formed by an 
intersection of said planar entrance and exit sides, and a point where said beam of 
light exits said exit side, determined by measurement and SnelTs-law calculation 
when said vessel contains the reference liquid for calibration, 
20 6^ 0 is a reference value of said angle with respect to the normal to the surface 

of the exit side at which said beam of light exits said exit side, determined when the 
vessel contains a reference liquid for calibration, 

D is distance from said point where said beam of light exits said exit side to 
said position sensor, determined when the vessel contains the reference liquid for 
25 calibration, 

W(n,) = is a value of said distance between said inner corner formed by the 
intersection of said planar entrance and exit sides, and a point where said beam of 
light exits said exit side, when said vessel contains the sample of the multi-component 
liquid mixture, defined according to a second equation: 
30 Wfa) = L cotfarc sin((sin OJ/n^] + 1 tan{arc sin[(n/ng)(l - (sine^n,)' 7 ']} 
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where: 

L - distance from said inner corner to a point where said beam of light 
passes from said entrance side into the given sample of the multi-component 
liquid mixture, determined by measurement and Snell's-law calculation, 
5 n, = refractive index of the sample of the multi-component liquid 

mixture which is to be determined, 

t = thickness of said entrance side and said exit side of said vessel, and 
n g = refractive index of said entrance side and said exit side of said 

vessel, 
10 and 

0 c (ni) is the value of 0 C , for the sample of the multi-component liquid mixture, 
defined according to a third equation: 

6e( n i) = arc sin {n, cos[arc sin ((sin 6/111)]}. 

4. The method of Claim 3, further comprising the additional step of performing 
1 5 calibration and measurement to determine W 0 , 8^ 0 D, said calibration and 
measurement including: 

causing a first reference liquid of known refractive index to be contained in 
said vessel; 

causing a relative rotation between said vessel and said beam position sensor 
20 such that, for the first reference liquid, said beam of light exits said exit side in a 

direction substantially perpendicular to said beam position sensor, whereby said beam 
of light impinges on the beam position sensor at the reference position; 
measuring W 0 and recording 0^ o ; 
emptying said vessel of the first reference liquid; 
25 causing a second reference liquid of known refractive index to be contained in 

said vessel, the refractive index of the second reference liquid being different than that 
of the first reference liquid; 

determining values of d,W and 9 e corresponding to the second liquid, the 
values of W and 0 C being respectively calculated from said second and third 
30 equations; and 
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using the values of d, W and 0 e in said first equation to determine D. 


5. The method of Claim 1, further comprising the additional step of recording the 
temperature of the sample of the multi-component liquid mixture, wherein step (b) 
comprises: 

5 (b- 1 ) comparing the refractive index of the given sample of the multi- 

component liquid mixture determined in step (a) to values of the predetermined data 
relating the different concentrations of the given component of the multi-component 
liquid mixture to corresponding values of the refractive index of the multi-component 
liquid mixture, for at least two temperatures near the recorded temperature of the 

1 0 sample of the multi-component liquid mixture; and 

(b-2) interpolating between the values of the predetermined data at the at least 
two temperatures near the recorded temperature of the sample of the multi-component 
liquid mixture so as to obtain an interpolated concentration of the sample of the multi- 
component liquid mixture. 

1 5 6. The method of Claim 1 , wherein the predetermined data relating different 
concentrations of the given component of the multi-component liquid mixture to 
corresponding values of the refractive index of the multi-component liquid mixture 
are double valued, such that individual values of the refractive index in at least a 
portion of a given range of the corresponding values of the refractive index 

20 correspond to two different values of concentration, said method further comprising 
the additional steps of: 

(c) storing previously measured values of concentration for the given sample of 
the multi-component liquid mixture taken at prior points in time, said previously 
measured values of concentration being indicative of the expected range of values of 

25 the concentration being determined; and 

(d) using the previous values of the concentration stored in step (c) to select which 
of the two different values of the concentration corresponding to the value of the 
refractive index determined in step (a) is correct. 
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7. A method for determining a change in concentration of a given component in a 
sample of a multi-component liquid mixture from an initial concentration of the given 
component corresponding to an initial refractive index of the sample at an initial 
temperature, said method comprising the steps of: 
5 (a) determining the refractive index of the sample of the multi-component liquid 
mixture by: 

(a-1) causing the sample of the liquid mixture to be contained in a vessel 
having a planar entrance side and a planar exit side, with a known angular relationship 
therebetween, said vessel being immersed in known surroundings; 
1 0 (a-2) causing a beam of light to impinge on said entrance side at an angle 6j 

with respect to a normal to an outer surface of said entrance side, to pass through the 
sample of the multi-component liquid mixture, and to then pass through said exit side, 
from which it exits at an angle 0 e with respect to a normal to an outer surface of said 
exit side; and 

1 5 (a-3) calculating the refractive index of the sample of the multi-component 

liquid mixture, based on said angles 6j and 8 C , by applying Snell's law at: 

an interface between the surroundings and said entrance side; 
an interface between said entrance side and the given sample of 
the multi-component liquid mixture; 
20 an interface between the sample of the multi-component liquid 

mixture and said exit side; and 

an interface between said exit side and the surroundings: 
(b) determining a change in refractive index, An, by subtracting from the initial 
refractive index the refractive index determined in step (a-3); 
25 (c) measuring the temperature of the sample of the multi-component liquid 
mixture; and 

(d) determining the change, AC, in the concentration of the given component in 
the sample of the multi-component liquid mixture from the initial concentration 
according to the approximate formula: 
30 AC s (An - (an/3T)AT)(an/aC)-\ 

where: 
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dn/ffT is partial derivative of index of refraction, n, with respect to 
temperature, determined from known data, 

AT is the difference between the temperature measured in step (c) and the 
initial temperature, and 
5 dn/dC is partial derivative of index of refraction, n, with respect to 

concentration, evaluated from known data in a region near the initial concentration. 

8. A method for determining the temperature of a sample of a liquid, the liquid 
being one of a pure liquid and a multi-component liquid mixture having a 
substantially constant concentration of its components, said method comprising the 
10 steps of: 

(a) determining the refractive index of the sample of the liquid by: 

(a-1) causing the given sample of the liquid to be contained in a vessel having 
a planar entrance side and a planar exit side, with a known angular relationship 
between said entrance side and said exit side, said vessel being immersed in known 
15 surroundings; 

(a-2) causing a beam of light to impinge on said entrance side at an angle 6; 
with respect to a normal to an outer surface of said entrance side, to pass through the 
sample of the liquid, and to then pass through said exit side, from which it exits at an 
angle 0 C with respect to a normal to an outer surface of said exit side; and 
20 ( a -3) calculating the refractive index of the given sample of the liquid, based 

on said angles 6 S and 8^ by applying Snell's law at: 

an interface between the surroundings and said entrance side; 
an interface between said entrance side and the sample of the 

liquid; 

25 an interface between the sample of the liquid and said exit side; 

and 

an interface between said exit side and the surroundings; 

and 

(b) comparing the refractive index of the given sample of the liquid determined in 
30 step (a) to predetermined data relating different temperatures of the liquid to 
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corresponding values of the refractive index of the liquid, so as to determine the 
temperature of the sample of the liquid corresponding to the refractive index of the 
sample determined in step (a). 

9. The method of Claim 8, wherein: 

step (a) further comprises the additional step of: 

(a-4) providing a beam position sensor which senses a position of said 
beam of light after said beam exits said exit side, and wherein step (a-3) uses 
said sensed position of the beam in calculating the refractive index of the 
sample of the multi-component liquid mixture. 

1 0. The method of Claim 9, wherein: 

said entrance and exit sides of said vessel have a substantially uniform 
thickness and a substantially uniform refractive index; and 

step (a-3) comprises determining the refractive index of the sample of the 
liquid mixture according to a first equation: 
d(n,) = [W 0 -W(n,)] cos 8^ + {D-tWo-W^)]} tan^n,)^] 
where: 

d(n,) is said position where said beam of light impinges on said beam position 
sensor, expressed as a distance relative to a reference beam position, 

W 0 is a reference value of distance between an inner corner formed by an 
intersection of said planar entrance and exit sides, and a point where said beam of 
light exits said exit side, determined by measurement and SnelFs-law calculation, 

0^0 is a reference value of said angle at which said beam of light exits said exit 
side, determined when the vessel contains a reference liquid for calibration, 

D is distance from said point where said beam of light exits said exit side to 
said beam position sensor, determined when the vessel contains the reference liquid 
for calibration, 

W(n0 = is a value of said distance between said inner corner formed by the 
intersection of said planar entrance and exit sides, and a point where said beam of 
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light exits said exit side, for the sample of the liquid, defined according to a second 
equation: 

W(n,) = L cot[arc sin((sin O^/n,)] + 1 tan{arc sin[(n/ng)(l - (sinO^/ii,)*]} 

where: 

5 L = distance from said inner corner to a point where said beam of light 

passes from said entrance side into the given sample of the liquid, determined 
by measurement and Snell Vlaw calculation, 

n } - refractive index of the sample of the liquid which is to be 
determined, 

10 t = thickness of said entrance side and said exit side of said vessel, and 

rig = refractive index of said entrance side and said exit side of said 

vessel, 

and 

6 c (n,) is the value of 0^ for the sample of the liquid, defined according to a 
15 third equation: 

8 e (n,) = arc sin {n, cos [arc sin ((sin 0/r^)]}. 

1 1 . The method of Claim 1 0, further comprising the additional step of performing 
calibration and measurement to determine W 0 , 0^ o and D, said calibration and 
measurement including: 
20 causing a first reference liquid of known refractive index to be contained in 

said vessel; 

causing a relative rotation between said vessel and said position sensor such 
that, for the first reference liquid, said beam of light exits said exit side in a direction 
substantially perpendicular to said position sensor, whereby said beam of light 
25 impinges on the beam position sensor at the reference position, whereby the value of d 
for the beam at said reference position may be taken as zero; 

measuring W 0 and recording 0 C 0 ; 

emptying said vessel of the first reference liquid; 
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causing a second reference liquid of known refractive index to be contained in 
said vessel, the refractive index of the second reference liquid being different from 
that of the first reference liquid; 

determining values of d,W and 6 C corresponding to the second liquid, the 
5 values of W and G e being respectively calculated from said second and third 
equations; and 

using the values of d, W and 0 e in said first equation to determine D. 

12. An apparatus for detennining the concentration of a given component in a 
sample of a multi-component liquid mixture, said apparatus comprising: 
1 0 (a) a vessel which contains said sample, said vessel having an entrance side and 
an exit side, with a known angular relationship between said entrance side and said 
exit side, said vessel being immersed in known surroundings; 

(b) a light source which is located so as to cause a beam of light to impinge on 
said entrance side at an angle Q t with respect to a normal to an outer surface of said 

1 5 entrance side, to pass through said sample, and to then pass through said exit side, 
from which it exits at an angle 6 C with respect to a normal to an outer surface of said 
exit side; 

(c) a sensor which determines said exit angle 0 C , thereby permitting the refractive 
index of said sample to be determined, based on said angles 0; and 8 e , by applying 

20 Snell's law at: 

an interface between the surroundings and said entrance side; 

an interface between said entrance side and the sample; 

an interface between the sample and said exit side; and 

an interface between said exit side and the surroundings; and 
25 (d) a computer programmed to compare a refractive index of said sample 
determined using said angles 0 S and 6 e with predetermined data relating different 
concentrations of the given component of the multi-component liquid mixture to 
corresponding values of the refractive index of the multi-component liquid mixture so 
as to determine the concentration of the given component in the sample of the multi- 
30 component liquid mixture.. 
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1 3 . The apparatus of Claim 1 2, wherein: 

said sensor comprises a beam position sensor which senses a relative position 
of said beam of light after said beam exits said exit side, said beam position sensor 
having a known geometrical relationship to said vessel such that said sensed beam 
5 position can be related to said exit angle 0 C . 

1 4. The apparatus of claim 12 wherein the computer stores the predetermined data 
relating different concentrations of the given component of the multi-component 
liquid mixture to corresponding values of the refractive index of the multi-component 
liquid mixture. 


10 15. The apparatus of Claim 12, wherein said entrance and exit sides of said vessel 
have a substantially uniform thickness and a substantially uniform refractive index, 
said apparatus and said computer receives a representation of said position of said 
beam of light and which determines the refractive index of the material according to a 
first equation: 

15 d(n,) - [W 0 -W(n>)] cos 0 cO + {D-[W 0 -W( ni )]} tan^n,)^] 

where: 

d(n,) is said relative position of said beam of light, measured with said position 
sensor, relative to an initial calibration value, 

is a reference value of distance between an inner corner formed by an 
20 intersection of said entrance and exit sides, and a point where said beam of light exits 
said exit side, determined during calibration, 

0^ o is a reference value of said angle at which said beam of light exits said exit 
side, determined during calibration, 

D is distance from said point where said beam of light exits said exit side to 
25 said position sensor, determined during calibration, 

W(n,) is a value of said distance between said inner corner formed by the 
intersection of said entrance and exit sides, and a point where said beam of coherent 
light exits said exit side, for the sample, during a given measurement, defined 
according to a second equation: 
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W(n,) = L cotfrrc sin((sin 0^)] + 1 tan{arc sin^n/n^O - (sin 

where: 

L = distance from said inner corner to a point where said beam of 
coherent light passes from said entrance side into the liquid, determined during 
calibration, 

n, = refractive index of sample which is to be determined, 
t = thickness of said entrance side and said exit side of said vessel, and 
ng = refractive index of said entrance side and said exit side of said 

vessel, 

and 

0 c (n,) is the value of 6 C , for the liquid, during said measurement, defined 
according to a third equation: 

0 c (n,) = arc sin {ni cos[arc sin ((sin QJn x )]}. 

15 16. The apparatus of Claim 12, wherein: 

said sensor comprises a photodiode beam position sensor and a position-to- 
voltage converter which produces an output voltage which is substantially linearly 
proportional to said relative position of said beam; and 
said computer is a digital computer; 
20 said apparatus further comprising: 

a first voltmeter which receives said output voltage from said position-to- 
voltage converter and which supplies a corresponding digital signal to said computer 
as said representation of said relative position of said beam of light; and 

a temperature sensor which is positioned to sense temperature of the sample of 
25 the liquid mixture and which is coupled to said computer so as to supply said 

computer with a digital signal representative of the temperature of the sample of the 
liquid mixture. 

1 7. The apparatus of Claim 1 6, further comprising a rotation stage which is 
interposed between said vessel and a mounting surface such that said relative position 
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of said beam of light can be adjusted, during calibration of said apparatus, to be 
substantially perpendicular to said beam position sensor and impinging on said beam 
position sensor at a reference position thereof. 

1 8. The apparatus of Claim 1 7, further comprising a scale interposed between said 
5 rotation stage and the mounting surface, whereby changes in mass concentration of 
the liquid may be determined. 


19. An apparatus for determining the temperature of a sample of liquid, said 
apparatus comprising: 

(a) a vessel which contains said sample, said vessel having an entrance and an 
1 0 exit side, with a known angular relationship, between said entrance side and said exist 

side, said vessel being immersed in known surroundings; 

(b) a light source which is located so as to cause a beam of light to impinge on 
said entrance side at an angle 0, with respect to a normal to an outer surface of said 
entrance side, to pass through said sample, and to pass through said exit side, from 

1 5 which it exits at an angle 8 C with respect to a normal to an outer surface of said exit 
side; 

(c) a sensor which determines said exit angle 6 e thereby permitting the 
refractive index of said sample to be determined, based on said angle 0 ; and 9 C , by 
applying SnelFs law which: 

20 an interface between the surrounding and said entrance side; 

an interface between said entrance side and the sample; 
an interface between the sample and said exit side; and 
an interface between said exit side and the surroundings; 

(d) a computer programmed to compare a refractive index of said sample 
25 determined using said angles Q { and 0 C with predetermined data relating different 

temperatures of the liquid to corresponding values of the refractive index of the liquid 
so as to determine the temperature of the sample of the liquid. 
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20. The apparatus of Claim 1 9 wherein the computer stores the predetermined data 
relating different temperatures of the liquid to corresponding values of the refractive 
index of the liquid. 

21 . The apparatus of claim 12, wherein said computer is programmed to perform 
5 the comparison of the determined refractive index to values of the predetermined data 

for at least two temperatures near the temperature of the liquid mixture and to 
interpolate between the values of the predetermined data at the at least two 
temperatures so as to correct the value of the concentration of the given component in 
the sample of the multi-component liquid mixture for the temperature of the sample of 
10 the liquid mixture. 

22. The apparatus of claim 12, wherein said computer stores previously measured 
values of concentration of the given component in the sample of the multi-component 
liquid mixture taken at prior points in time and said computer is programmed to 
employ the previously measured values to indicate an expected range of 

1 5 concentrations of the given component to permit the selection of a correct value of 
concentration when concentration is double-valued with respect to index of refraction. 


23. The apparatus of Claim 12, further comprising: 

a temperature sensor which is positioned to sense temperature of the sample of 
the multi-component liquid mixture and which is coupled to said computer so as to 
20 supply said computer with a signal representative of the temperature of the multi- 
component liquid; 
wherein: 

said computer is programmed to determine the refractive index of a given 
sample of the multi-component liquid mixture, based on said angles Q { and 6 C , by 
25 applying Snell's law at: 

an interface between the surroundings and said entrance side; 
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an interface between said entrance side and the sample of the 
multi-component liquid mixture; 

an interface between the sample of the multi-component liquid 
mixture and said exit side; and 
5 an interface between said exit side and the surroundings; 

and; 

said computer is also programmed to determine a change in refractive index of 
the multi-component liquid mixture, An, by subtracting, from the refractive index 
determined by said computer using SnelFs law, an initial refractive index of the 

1 0 sample of the multi-component liquid mixture at an initial concentration of the given 
component in the sample and at an initial temperature of the sample measured with 
said temperature sensor, to receive said signal which represents the temperature of the 
sample of the multi-component liquid mixture under test, and to determine the 
change, AC, in the concentration of the given sample of the multi-component liquid 

1 5 mixture from the initial concentration according to the approximate formula: 
AC s (An - (an/aT)AT)(3n/aC)-\ 

where: 

dn/5T is partial derivative of index of refraction, n, with respect to 
temperature, determined from known data, 
20 AT is the difference between the temperature of the sample of the multi- 

component liquid mixture under test and the initial temperature, and 

dnfdC is partial derivative of the index of refraction, n, with respect to 
concentration, evaluated from known data in a region near the initial concentration. 
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